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Abstract. We present high-speed long term photometric ob- 
servations from 1992 to 1996 (UBVRI), as well as recently 
obtained (1997) spectrophotometric and high resolution spec- 
troscopic studies of the bright galactic supersoft X-ray source 
RX J0019.8 +2156. Our photometry reveals a highly variable 
object. The shortest observed quasi-periodic variations are 
humps with a period of «1.8 h. The timings of the main mini- 
mum are not exactly regular and occur on average around phase 
4> — 0.0. With our new data set we calculated an update of the 
orbital period. 

Our simultaneous spectroscopic and spectrophotometric 
studies from 1997 give an even more detailed insight: blue 
and red shifted satellite lines at He II A 4686 and at the strong 
Balmer lines, which are interpreted as high velocity outflows, 
can clearly be detected. We also present for the first time an 
analysis of this supersoft source by means of Doppler tomog- 
raphy. The lines of He II ( AA 4542, 4686), U a and Up show no 
emission from an accretion disk. These emission lines originate 
in regions within the binary system with very low velocities. 
The spatial origin of the emitting material is not quite clear and 
will be investigated in a follow-up paper. 

Key words: stars individual: QR And, - accretion disks - bi- 
naries: close - binaries: eclipsing - X-rays: stars 



1. Introduction 

Supersoft X-Ray sources (hereafter called SSS) were discov- 
ered by Einstein, established by Rosat (Trumperet al. 1991) 
and are now commonly accepted to be close binary systems 
with a white dwarf primary. Little is known about the nature of 
the companion star. The accretion rate in these systems is rel- 
atively high (in the order of 10 MQy -1 ) so that the accreted 
hydrogen burns steadily on the surface of the white dwarf (van 
den Heuvel et al. 1992). Several SSS are now known, most of 
them in the LMC and in M3 1 . A detailed list can be found in 
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Greiner (1996). Only two SSS are known in our Galaxy; one of 
them is RXJ0019. 8+2156 (hereafter called RXJ0019), which 
was discovered by Beuermann et al. (1995). 

SSS are characterized by their high X-ray luminosity from 
10 36 erg s _1 to 10 38 erg s _1 and by extremely soft X-ray spec- 
tra with blackbody temperatures between 20-60 eV. Most of 
the optical luminosity probably originates in the bright accre- 
tion disk: in a radially extended and vertically elevated accre- 
tion disk rim a substantial fraction of the soft X-ray flux from 
the white dwarf is reprocessed to optical and UV-light. This 
model successfully accounts for features observed in the op- 
tical lightcurves of some SSS such as CAL 87, RXJ0019 and 
RXJ0513.9-6951 (Meyer-Hofmeister et al. 1997). Most of the 
optical lightcurves presented in this paper have also been used 
by Meyer-Hofmeister et al. (1998) for their analysis . They ar- 
gue that the short time variability is caused by changes in the 
height of the accretion disk rim. An elevated rim is therefore 
the premise for explaining the luminosity in optical bands. This 
picture of a bright accretion disk is also invoked by other au- 
thors (Popham & Di Stefano 1996, Matsumoto & Fukue 1998). 

RXJ0019 is a unique object. It is the brightest and nearest 
among the known SSS. With its position well outside the galac- 
tic plane and its visual brightness of my w 12.5™ RXJ0019 
offers an outstanding opportunity to investigate a SSS very 
closely with indirect imaging techniques such as Doppler to- 
mography. Here we report on such an approach to investigate 
RXJ0019. 

In Sect. 2 we refer to our observations of RXJ0019 and 
the data reduction procedures used. Section 3 deals with the 
results and analysis of our data. Finally, in Sect. 4 a summary 
and discussion of this paper are given. 

2. Observations and data 

RXJ0019 has been monitored regularly at the Wendelstein Ob- 
servatory from 1992 to 1995 (Will & Barwig 1996). Follow-up 
observations were made in autumn 1996. We were able to ob- 
serve this object on 15 nights between October and December 
1996. 

A special approach to investigate RXJ0019 was made in 
October 1997. We performed simultaneous spectroscopy and 
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Table 1. Journal of photometric observations with MCCP. 
HJD= truncated Heliocentric Julian Date -2440000.0, 
h = hours, IT = integration time 



Date 


Start [HJD] 


Stop [HJD] 


Duration [h] 


IT [s] 


21/09/1992 


08887.296 


08887.668 


8.93 


2 


22/10/1995 


10013.225 


10013.648 


10.16 


2 


23/10/1995 


10014.249 


10014.637 


9.31 


2 



2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
1 
2 
2 
2 
2 



spectrophotometry at the 3.5 m and 2.2 m telescope at Calar 
Alto, respectively. Detailed information about the individual 
observing campaigns is given in Tables to | 

2.1. Photometry 

High-speed photometric observations of RXJ0019 were per- 
formed using the multichannel multicolor photometer MCCP 
(Barwig et al. 1987) at the 80 cm telescope at the Wendelstein 
Observatory in the Bavarian Alps. With this photometer we 
can monitor the object, a nearby comparison star and the sky 
background simultaneously in UBVRI. The MCCP allows a 
nearly complete elimination of atmospheric transparency vari- 
ations and extinction effects (Barwig et al. 1987), by subtract- 
ing the sky background of each color channel and dividing the 
object by the comparison star measurement afterwards. There- 
fore, photometric measurements were possible even under non- 
photometric conditions. During the 1996 run we were able to 
add another 88.7 hours of observation time for RXJ0019 to our 
existing dataset. 

As comparison star we used the star at RA — 19 m 31!2, 
DEC = +21°53'58'/9 (for epoch 2000.0). The intensities of 
the light-curves presented in Fig. [3] are calculated relative to 
this star. The error in the relative count rate is ±0.012. In order 
to normalize the individual channels, calibration measurements 
were performed during photometric conditions once or twice a 
night. The constancy of the calibration coefficients is indica- 
tive of the high stability of our detectors. All observations were 
made with an integration time of 2 sec. except on the night of 
Dec 10,1996 where we used an integration time of 1 sec. A de- 
tailed journal of the observations is given in Table [l] 



Table 2. Journal of spectrophotometric observations 
with MEKASPEK. HJD = truncated Heliocentric Julian Date 
-2450000.0, h = hours, IT = integration time 



Date Start [HJD] 


Stop [HJD] 


Duration [h] 


IT[s] 


27/10/1997 


749.418 


749.628 


4.23 


2 


30/10/1997 


752.295 


752.623 


7.74 


2 


Table 3. Journal of spectroscopic 


observations. 


B=blue 


spectral region. 


, R=red spectral region, UT=universal time, 


IT=integration time 








Date 


l ) UC-l, Ll Cl 


Phase range 


Start [UT] 


IT[s] 




I> 


0.12-0.15 


22 h 22 m 30 s 


1800 




B R 


0.16-0.19 


23 h ll m 55 s 


1800 




B R 


0.20-0.23 


23 ,i 48 m 25 s 


1800 


27-28/10/1997 


B R 


0.24-0.27 


00 ,i 23 m 42 s 


1800 




B R 


0.28-0.31 


01 h 00 m 01 s 


1800 




B R 


0.33-0.37 


0l' l 50 m 13 s 


1800 




B 


0.37-0.40 


02 h 25 m 04 s 


1800 


28-29/10/1997 


B R 


0.69-0.72 


22 h 56 m 57 s 


1800 




B R 


0.44-0.47 


18 h 54 m 08 s 


1800 




B R 


0.48-0.51 


19 h 28 m 49 s 


1800 




B R 


0.51-0.54 


20 ,1 03 m 14 s 


1800 




B R 


0.55-0.58 


20 ,1 37 m 39 s 


1800 




B R 


0.59-0.61 


21 h 21 m 18 s 


690 




B R 


0.69-0.72 


22 h 52 m 35 s 


1800 


30-31/10/1997 


B R 


0.73-0.76 


23 ,t 26 m 55 s 


1800 




B R 


0.76-0.79 


00' l 01 m 13 s 


1800 




B R 


0.80-0.83 


00 ,l 35 m 47 s 


1800 




B R 


0.83-0.86 


01 h 10 m 30 s 


1800 




B R 


0.87-0.90 


01 h 46 m 23 s 


1800 




B R 


0.91-0.94 


02 fc 20 m 45 s 


1800 




B R 


0.95-0.98 


02 h 55 m 20 s 


1800 



2.2. Spectrophotometry 

We observed RXJ0019 on October 27 and 30, 1997 with 
MEKASPEK attached to the 2.2 m telescope at the Calar 
Alto Observatory. MEKASPEK is a four channel fiber-optic 
spectrophotometer, developed at the Universitats-Sternwarte 
Miinchen. With MEKASPEK we can also perform simultane- 
ous measurements of the object, a nearby comparison star and 
the sky background within the spectral range of 3700...9000A 
at a spectral resolution of A/ A A w 50. The photon-counting 
two-dimensional detector (MEPSICRON) has a time resolution 
of up to 5 ms. 

Furthermore, MEKASPEK performs a correct treatment 
of atmospheric effects and allows an accurate transformation 
to any broadband photometric system. For more details see 
Mantel et al. (1993) and Mantel & Barwig (1993). Again, at- 
mospheric effects are eliminated using the standard reduction 
method (Barwig et al. 1987). 

We used an integration time of 2 sec. Each night two cal- 
ibration measurements were performed to get the calibration 
coefficients for the normalization of the color channels. We 
also made flat-field measurements and a measurement of a 
HgArRb-lamp spectrum for wavelength calibration at the end 



11/10/1996 10368.350 

12/10/1996 10369.261 

13/10/1996 10370.442 

18/10/1996 10375.321 

22/10/1996 10379.482 

23/10/1996 10380.237 

02/11/1996 10390.401 

03/11/1996 10391.224 

04/11/1996 10392.240 

09/11/1996 10397.218 

10/12/1996 10428.394 

11/12/1996 10429.248 

12/12/1996 10430.250 

15/12/1996 10433.210 

16/12/1996 10434.276 



10368.664 


7.537 


10369.655 


9.442 


10370.654 


5.111 


10375.552 


5.530 


10379.637 


3.704 


10380.476 


5.746 


10390.557 


3.630 


10391.614 


9.361 


10392.599 


8.597 


10397.507 


6.943 


10428.499 


2.515 


10429.417 


4.406 


10430.445 


4.697 


10433.480 


6.478 


10434.486 


5.051 
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Table 4. Magnitudes m max , m m i n and orbital modulation 
Am ofRXJ0019.8on Oct 12,1996. 



Minima 1 996 



Band 






Am 


U 


11.37 


11.79 


0.42 


B 


12.29 


12.71 


0.42 


V 


12.33 


12.65 


0.42 



of each night. Due to unfortunate weather conditions only two 
of four program nights could be used. Therefore, we could not 
cover a complete orbital phase of RXJ0019. The details of the 
MEKASPEK observations are listed in Table §. 

2.3. Spectroscopy 

Simultaneously with the spectrophotometric measurements 
we obtained high-resolution spectroscopic data with the 
Cassegrain double beam spectrograph (TWIN) attached to the 
3.5 m telescope at Calar Alto. The blue and red channel were 
equipped with the low-noise CCDs (SITe#6a and SITe#4d) 
with a pixel size of 15/im and a CCD size of 800x2000 
pixel. We used a slit width of 1'.'5. For the blue channel we 
chose grating T07 with a spectral resolution of 0.81 A/pixel 
and a spectral range of 3300...5000A (observations performed 
in second order) and for the red channel grating T04 with 
a spectral resolution of 1.08 A/pixel in a spectral range be- 
tween 5500...9000A. The object was trailed along the slit in 
order to get phase-resolved, high resolution spectra. Expo- 
sure times ranged between 690 and 1800s and the mean trail 
velocity was set to 320"/h. The trailed spectra were binned 
in the direction of the slit by a factor of 2. Helium- Argon 
wavelength calibration spectra were taken approximately ev- 
ery hour. During the first observing night no calibration spec- 
tra in the blue spectral range could be obtained due to tech- 
nical problems. Flatfield measurements were taken at the end 
of the last observing night. Further information about the ob- 
servation is given in Table [}| The data reduction included 
bias-subtraction, flatfield-correction, sky-subtraction, cosmic- 
ray elimination and wavelength-calibration as described by 
Home (1986). 

After that, the continuum of the spectroscopic data set was 
calibrated with our simultaneously obtained spectrophotomet- 
ric data. By applying this method we were able to correct for in- 
tensity variations of the emission lines which might have been 
caused by poor weather conditions (e.g. clouds) or observa- 
tional difficulties (e.g. variable vignetting of the seeing disk by 
the slit). Further information is given in Simic et al. (1998). In 
a last step the trailed spectra were phase-folded into 125 phase 
bins using the updated eclipse ephemeris as given in Sect. 3.1. 

3. Analysis and results 

3.1. The lightcurves 

All observed MCCP and MEKASPEK lightcurves are pre- 
sented in Fig. B. Only B-band lightcurves are plotted since 
the very blue continuum of RXJ0019 in combination with the 




Fig. 1. Minima observed in 1996. The minima occur statisti- 
cally around phase cf> — 0.0. The depth and the shape of the 
minima varies. All lightcurves were approximated by splines. 



Colour variations on 09.1 1 ,9& 



Fig. 2. Color Variations on Nov 9, 1996. Symmetrically to 
phase 0.0 an increase in the R — I flux can be observed with 
an amplitude of ss 0.15™. The zero point of the vertical scale is 
arbitrary. 

quantum efficiency of the detectors and the color temperature 
of the comparison provides the best S/N value in B. All obser- 
vations are plotted individually to emphasize the high variabil- 
ity of this object. A plot including all observed lightcurves su- 
perposed according to the orbital period can be found in Meyer- 
Hofmeister et al. (1998) in their Fig. 1. Also compare with Fig. 
2 in Will & Barwig (1996). The following features can be ob- 
served: 

1.) RXJ0019 shows a deep, broad minimum (Am = 0.42™ 
in UBV on Oct. 12, 1996, Table g), lasting about 0.55 of a com- 
plete orbital phase. As can be seen in Fig. |l]the depth and shape 
of the minima varies and the timings of the minima are not reg- 
ular. This explains the systematic difference of the orbital pe- 
riod of Will & Barwig (1996), which lies out of the error bars 
of other periods given in the literature (e.g. Greiner & Wenzel 
1995 or Matsumoto 1996): our dataset was too short to elim- 
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Fig. 3. MCCP and MEKASPEK B-band lightcurves between September 1992 to October 1997. Three lightcurves from 1992 
and 1995 have also been included (Will & Barwig, 1996). Note that in 1995 the brightness of RXJ0019 is slightly lower. All 
lightcurves differ significantly from each other. Also note the quasi-periodic 1.8h humps on Sept 21,1992, Oct 11,1996, Oct 
13,1996. Increased scatter in the lightcurves is due to very poor observing conditions. 



Table 5. Timings of primary minima (HJD) 



date 


cycle (E) 


minima 


mean 






HJD-2450000 [d] 


error [d] 


12/10/1996 


22061 


369.590 


0.00099 


22/10/1996 


22076 


379.504 


0.00175 


03/11/1996 


22094 


391.385 


0.00774 


09/11/1996 


22103 


397.333 


0.00065 


12/12/1996 


22153 


430.361 


0.00130 


16/12/1996 


22159 


434.323 


0.00036 



inate the statistical effects. With the epoch of the photometric 
minimum given by Greiner & Wenzel ( 1 995) and the very accu- 
rate determination of the minima observed at the Wendelstein 
observatory since 1992 we derive the following orbital period: 

HJD = 2435799.247 + 0.6604573(40) x E 



The timings of the minima used for this calculation can be 
found in Table | and in Will & B arwig (1996) in their Table 
2. 

2. ) A secondary minimum appears around phase <j> 0.5. 
The amplitude of the secondary minimum is Am w 0.1 m . The 
following maximum at <j> « 0.65 is slightly fainter than the 
first maximum. Sometimes the second maximum is even com- 
pletely smeared out and after a flat portion of the lightcurve 
the ingress to the main minimum begins (e.g. Oct 22, 1995 and 
Dec 15, 1996 in Fig. |). 

3. ) We observe humps and step-like features during ingress 
and egress of the main minimum. A period analysis of these 
features, when apparent in the lightcurves, reveals a period of 
approximately 1.8 h (e.g. Oct 11,1996 and Oct 13,1996). In 
contrast to the lightcurves showing steps and humps there are 
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also lightcurves with a rather flat appearance (e.g. Nov 09, 1996 
or Dec 12,1996). We assume that RXJ0019 changes between 
two states: an "excited optical state" where humps and steps 
appear quasi-periodically with a period of 1.8 h and a "quiet 
optical state", where the lightcurves show no additional fea- 
tures. RXJ0019 can change from one state into the other from 
night to night. A probable explanation might be a short-term 
variation in the mass accretion rate (Meyer-Hofmeister et al. 
1998). 

4.) Matsumoto (1996) did not observe any color variations 
in his lightcurves. We also cannot find any color variations in 
U - B, B - V and V ~ R. But in R - / we detect a weak flux 
increase symmetrically to <fi — 0.0 with an amplitude of Am 
k 0.15 m . This variation in R — I can only be observed when 
the main minimum appears at midnight hours. In the evening 
and in the morning the natural reddening of the sky prevents an 
observation of this color variation. Therefore, we observed this 
color variation only on Oct 12, Nov 3 and Nov 9, 1996 (Fig. 
||). In every night the width and the amplitude of this variation 
are similar. All other portions of the color curves from phase 
<j> = 0.1 to 0.9 are flat and show no variations. 




3700 3800 3900 4000 4100 

Wovelenght (A) 



Fig. 5. Phase-resolved, uncalibrated spectra showing the vari- 
ability of the emission lines between 3700 A and 4100 A. The 
first spectrum from = 0. 17 — 0.21 is dominated by P-Cygni 
absorptions at the Balmer lines up to H12, which cannot be seen 
at later phases. Marks at the bottom indicate the Balmer series 
up to H i3 . 



3.2. The mean orbital spectra 

High-resolution spectra of RXJ0019 have recently been ob- 
tained by Tomov et al. (1998) and Becker et al. (1998). 
We also present high resolution mean orbital spectra (in Fig. 
0) in the blue (3700...5000A) and in the red spectral range 
(6300.. .8300A). Similar to other SSS (e.g.Cramptonetal. 1996 
and Southwell et al. 1996) the very blue spectrum is dominated 
by He II and Balmer lines. The Balmer series can be traced up 
to H12. All transitions of the Hell (n,3),(n,4) and (n,5) series 
can be observed when not blended by Balmer lines. No He I 
can be detected indicating a high ionisation state. 

Additionally, we see some higher ionization emission fea- 
tures such as O VI A 381 1. The C III-N III A A 4640 -4660 emis- 
sion complex contributes to the blue wing of the He II A 4686. 

On the blue wings of the Balmer lines P-Cygni profiles are 
visible in all observed lines. The Hg to H12 transitions can 
only be detected by their corresponding P-Cygni absorption 
between = 0.12 to 0.30 (Fig. We do not observe these 
absorptions after <p = 0.30. At the other Balmer lines the P- 
Cygni absorption shows an orbital modulation that almost dis- 
appears between phase = 0.7 and 0.9 (Fig. ||), indicating that 
there might be a directed fast wind in this system. The velocity 
of the wind amounts to Rs590kms~ 1 with the blue absorption 
wing extending to w900kms _1 . The Hell emission lines are 
not truncated by P-Cygni profiles. 

Two emission features, which are clearly above the S/N ra- 
tio at AA 4500, 4930, could not be identified with any reason- 
able ion species. 

The Balmer lines H Q , Hp and H 7 have broad bases and 
show satellite lines symmetrically around the main emission. 
All bright Balmer lines are blended by He II emission lines, 
in particular the two Helium lines AA 6527, 6560 significantly 
contribute to H a and Hell A 4859 to H^. Therefore, substruc- 
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Fig. 6. Phase-resolved spectra showing the variability of the 
emission lines between 4500 A and 4900 A. The P-Cygni ab- 
sorption at Hp almost disappears between <fi = 0.7 — 0.9. Also 
note the change of the base of He II A 4686 on the red wing of 
the emission line. 



tures in these lines are blurred by the Helium emission, whereas 
in H 7 , H,5 and H c similar complex substructures of different 
components can be detected. For illustration the trailed H 7 
spectrum is included in Fig. [k]. 

In Hell A 4686 we only observe a broad base. The adja- 
cent satellite components are not clearly separated from the 
main emission and are asymmetric (Fig. |] and|6|). Tomov et al. 
(1998) and Becker at al. (1998) argue that these satellite lines 
are the spectral signatures of jets originating near the white 
dwarf. These collimated high velocity outflows have recently 
also been reported by Quaintrell & Fender (1998) from their 
infrared spectroscopy. 
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Fig. 4. Mean orbital spectrum in the blue spectral range range from 4000A to 5000A and in the red spectral range from 6400A to 
7250A. All Balmer lines show P-Cygni absorption on their blue wings. The base of He II A 4686 is very broad and asymmetric 
but the jet lines are not separated from the main emission. The jet lines (labeled Si and S2) at H Q and H^ are clearly visible. H Q 
is considerably blended by two He II emission lines. 

Transient jets seem to be a common feature among the 
SSS. Jets have also been detected in RXJ05 13.9-6951 (e.g. 
Southwell et al. 1996) and recently in RX J0925 .7-4758 (Motch 
1998). The observed projected high outflow velocities of these 
sources, which are seen at very low inclination, indicate an ori- 
gin of the jet near the white dwarf. Therefore, the jet velocity 
is close to the escape velocity of the central object (see Livio 
1998) and can be used for estimating the inclination of the sys- 
tem. 

The projected velocities of the jet lines (measured at max- 
imum of their emission profiles) of RX J0019.8 are rather low. 
We measure 920 kms -1 for Si and 805 kms -1 for S2- The 
asymmetry can be explained by the P-Cygni absorption, which 
alters the appearance of the jet emission on the blue side sig- 
nificantly. The low velocities may be due to a high system in- 
clination (Tomov et al. 1998) if the jet originates near the white 
dwarf. A medium to high inclination is also consistent with the 
deep eclipse lightcurves from our photometry. 

We note, that no spectral features of a secondary star are 
detected in our high resolution spectra. 




3.3. Radial velocities 

The usual procedure of measuring radial velocities in complex 
emission lines is fitting Gaussians to their line profiles. Measur- 
ing radial velocities of the emission lines of the RXJ0019 spec- 
tra is difficult. First, the Balmer emission lines are truncated 
by the phase dependent P-Cygni profiles and blended by He II 
emission lines which will alter the system velocity and semi 
amplitude considerably. Second, the asymmetric profile of the 
satellite lines will cause systematic errors in the velocity deter- 
minations. Furthermore, the measured velocities are only of the 
order of the spectral resolution of the spectrograph. Therefore, 
we restricted our investigation to the strongest emission lines 
H a and He II. 

For the pixel rows of the central parts of our trailed spec- 
tra we fitted Gaussians to the emission line profile. In order to 
suppress the influence of the noise and the satellite lines, the 



0.4 0.6 

Phase 

Fig. 7. Radial velocity curve determined by fitting Gaussians to 
the trailed spectra of H Q . The upper panel shows the residuals, 
the lower the distribution of the data. The solid line is the sine- 
fit to the data. For comparison the radial velocities from Tomov 
et al. (1998) are also plotted as asterisks (different symbols of 
our measurements refer to different trailed spectra). 

central region of the Gaussians was given a bigger weight than 
the adjacent spectral range (Bevington & Robinson 1992). Af- 
ter that a sine-fit was applied to the fitted radial velocity V and 
the observed semi- amplitude K b s , using the following rela- 
tion: 



V = 7 + K obs sin(0 - 0o ) 



(1) 



where 7 is the system velocity, K obs is the radial semi- 
amplitude of the white dwarf (assuming these emission lines 
originate near the primary component), and <fi and </> indicate 
the phase and the phase offset of the binary system, respec- 
tively. 

As we do not have blue calibration spectra from the first 
observing night our velocity determinations of He II were re- 
stricted to the spectra obtained in the last night. 
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Fig. 8. Radial velocity curve determined by fitting Gaussians to 
the trailed spectra of He II. The upper panel shows the residu- 
als, the lower the distribution of the data. The solid line is the 
sine-fit to the data. As the CCD is slightly warped we measure a 
systematic shift of the radial velocities within a single spectrum 
(see text) to higher velocities. 




0.0 0.5 1.0 1.5 2.0 



Fig. 9. Primary mass of RXJ0019 for various inclinations, i, 
plotted from the mass function f(M) = 0.0274AfQ. For any 
reasonable inclination it is clear that the companion must be a 
low-mass star. 

There is a systematic shift in our velocity measurements: 
measurements from a trailed spectrum show a continuous drift 
to higher velocities. The reason for this is that the CCD chips 
are slightly warped. We account for this effect by fitting sky 
lines in the spectra from which we calculate the distortion. But 
the effect could not be eliminated in the blue spectral range, as 
there were only two weak sky lines which could not be used for 
a reasonable fit. Therefore, no correction was made in the blue 
spectral range. It should be noted that in spite of this systematic 
error our residuals are still low (Figs. ^ and 

The sine-fits to the radial velocity curves of H Q and He II 
are given in Figs. and §1 respectively. The deduced velocities 
are listed in Table |[ It is obvious, that there is a substantial 



Table 6. Orbital parameters for RXJ0019. K b s is the ampli- 
tude of the radial velocity and 7 denotes the system velocity in 
kms -1 . <fio denotes the phase offset between photometric and 
spectroscopic phase (see Text), a designates the errors of the 
mentioned parameters. The system velocities are not consistent 
within their error bars. See text for an explanation. 
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°K obs 


7 


Cry 


00 




Hell 


13.7 


2.6 


-131.7 


25.0 


0.11 


0.03 


H Q 


33.0 


1.8 


-53.0 


12.0 


0.17 


0.01 



difference between the radial velocities in H Q and He II: the 
semi-amplitude in He II is considerably higher. The difference 
in the radial velocity between the two lines can be explained 
by the blurred H Q emission (see Sect. 3.2). Our result for the 
radial velocity of Hell A4686 is 73.7 ± 2.6 kms -1 , which is 
in accordance with other values given in the literature (Becker 
et al. 1998, Beuermann et al. 1995). 

We also observe a considerable difference in the system ve- 
locity between H Q and He II in our data. This effect is caused 
by the P-Cygni profiles: the blue wing of the H Q emission is 
truncated by the P-Cygni absorption and shifts the maximum of 
emission to slightly longer wavelengths. Thus, the system ve- 
locity in H Q differs notably from the velocity measured in He II. 
Compared with other values published in literature (Becker et 
al. 1998) our derived system velocities are much higher, es- 
pecially for He II. For a reliable determination of the system 
velocity a dataset covering only a single orbital period is prob- 
ably too short to eliminate statistical effects, which are caused 
by the varying influence of the jet emission. For comparison 
we show in Fig. ^ the velocities given in Tomov et al. (1998) 
in their Table 1 . It is obvious that their system velocity for H Q 
is even higher. But due to the systematic shift in our distorted 
spectra our measured system velocity should not be given too 
much weight. 

Photometric minimum and spectroscopic phase zero do not 
coincide. The photometric minimum occurs 0.11 and 0.17 in 
phase before spectroscopic phase zero for He II and H a , re- 
spectively. A similar shift for He II was already observed by 
Beuermann et al. (1995). 

For further investigations we adopt our results of the semi- 
amplitude for the rest of this paper. The system velocities for 
the Doppler tomography are taken from Becker et al. (1998). 

As the jet lines of H Q mimic the orbital velocity modulation 
of He II (Fig. |l^), we assumed that the radial velocity modula- 
tion of He II is related to the primary and therefore represents 
the motion of the white dwarf. Knowing the orbital period and 
amplitude of the radial velocity variations of the mass accretor 
we can deduce the mass function of the secondary star. For the 
calculation we use the data from the He II emission line. The 
mass function is determined by 

P orb Kf M?sin 3 i , . 

MM,) = = = O.O274(28)M (2) 
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where P or b is the orbital period, K\ is the radial velocity of the 
primary, Mi and M2 denote the mass of the white dwarf and 
the secondary, respectively, and i is the inclination. 

The mass function is obviously very low. If the emission 
lines of He II are indeed related to the compact object we can 
deduce the mass of the companion. In Fig. ^ we plot our de- 
rived mass function for various inclinations i. Because of the 
arguments given in Sect. 3.2, we assume that the inclination 
is reasonably high. We conservatively estimate that the incli- 
nation ranges between 50° < i < 90°. If the compact object 
is a white dwarf with 0.6 < M% < 1.2Mq, then we obtain 
0.3 < M2 < 0.5Mq for the mass of the companion. The com- 
panion star is almost certainly a low-mass object. This could 
explain why we do not see any spectral features of the compan- 
ion in our high resolution spectra. As low mass functions are 
also known from other SSS (e.g. Crampton et al. 1996, Hutch- 
ings et al. 1998), a low mass companion star may be common 
among this class of objects. A possible evolutionary scenario 
for low-mass secondaries in SSS is proposed by vanTeeseling 
& King (1998). 

3.4. Doppler tomography 

Doppler tomography is an useful tool to extract further infor- 
mation on the emission line origin from trailed spectra. This 
indirect imaging technique, which was developed by Marsh & 
Home (1988), uses the velocity of emission lines at each phase 
to create a two-dimensional intensity image in velocity space 
coordinates (Vx, Vy)- The Doppler map can be interpreted as 
a projection of emitting regions in accreting binary systems 
onto the plane perpendicular to the observer's line of sight. The 
Doppler map is a function of the velocity (Vx, Vy), where the 
X-axis points from the white dwarf to the secondary and the 
F-axis points in the direction of the secondary's motion. 

An image pixel with given velocity coordinates (Vx, Vy) 
produces an S-wave with the radial velocity 

V = 7 - V x cos(2tt</>) + Vy sin(27T0) (3) 

where 7 and <fi denote the system velocity and the phase, re- 
spectively. 

To accomplish this, a linear tomography algorithm, the 
Fourier-filtered back-projection (FFBP) is used, which is de- 
scribed in detail by Home (1991). The resulting Doppler map 
(or tomogram) is displayed as a grey-scale image. To assist in 
interpreting Doppler maps we also mark the position of the sec- 
ondary star and the ballistic trajectory of the gas stream. As- 
suming a Keplerian velocity field, emission originating in the 
inner parts of the disk has a larger velocity and thus appears in 
the outer regions of the map. Therefore, an image in such a rep- 
resentation is turned inside-out. For the tomogram analysis we 
restricted our investigation to the two Balmer lines H a and H^ 
and the emission lines Hell (A 4686) and Hell (A 4542) . No 
reasonable Doppler maps could be produced for other emission 
lines because of their poor S/N ratios. 

Before computing the tomograms we subtracted the under- 
lying continuum from the individual emission lines, since their 



line flux is the quantity needed to produce these maps. This 
was done by subtracting from each pixel row of the trailed 
spectra the corresponding median of the intensity. We achieved 
a FWHM resolution of 150 kms -1 in the central regions of 
the tomograms, whereas at higher velocities (« 800 kms -1 ) the 
resolution of the maps suffers a considerable degradation. Lin- 
ear structures in the high velocity regions (Fig. [fl]) of the maps 
are sampling artifacts (aliasing streaks) which are not taken into 
account in the further interpretation. For a detailed discussion 
about sampling artifacts see e.g. Robinson et al. (1993). 

In Fig. |l0] we present the trailed spectra of the emission 
lines in V, ^-coordinates; the Doppler tomograms are shown in 
Fig. |ll]. The schematic overlays in the Doppler maps represent 
the Roche-lobe of the companion star and the gas stream. They 
are a function of K\ and K%. From our discussion in Sect. 3.3 
we have taken 1 Mq for the white dwarf and a 0.5A/q compan- 
ion star. The center of mass and the location of the white dwarf 
are respectively marked by a cross and a point below the Roche 
lobe. The ballistic stream is represented by an arc originating 
from the secondary's Roche lobe at the inner Lagrangian point. 
This arc is marked every Q.1Rl 1 (open circles, Rj J1 is the dis- 
tance from the center of the primary to the inner Lagrangian 
point) as it accelerates towards the compact object. 

In the Balmer emission lines the blue side is truncated by 
the P-Cygni absorption. No S-wave structures can be detected 
at H Q and Hp. The jet lines of H a clearly show an orbital mo- 
tion which is in phase with the motion of He II. The main emis- 
sion component of H Q reveals almost no orbital motion (com- 
pare Table ||). Therefore, we assume that the jet and the main 
emission do not have the same spatial origin. 

The trailed spectra of both He II lines show more details: in 
each spectrogram a weak S-wave component within the main 
emission can be detected. 

The Doppler tomograms (Fig. [FJ]) of our data immediately 
reveal that the emission line distribution does not resemble a 
typical accretion disk. The usual shape of an accretion disk 
would appear as a dark extended ring in the inverse grey-scaled 
Doppler images, e.g. compare with Doppler maps of IP Peg ob- 
tained with the same technique in Wolf et al. (1998). Our ob- 
served projected velocities are smaller than 350 kms -1 . For a 
binary system with parameters mentioned above and a disk ra- 
dius of r = 0.8Rl = 10 10 95 cm (Rl is the volume radius of 
the Roche lobe of the white dwarf) one finds a Kepler veloc- 
ity of 384 kms -1 at the rim of the disk. As the Kepler velocity 
increases rapidly within the accretion disk, the observed emis- 
sion cannot originate within the disk if the inclination is high. 
For medium inclinations an origin only at the rim of the disk 
is possible. This conclusion is valid for a huge variety of sys- 
tem parameters (Fig. |l2|). Low inclinations are excluded (see 
discussion in Sect. 3.2). 

The dim ring in the map of H a at a velocity of rs 950 
kms -1 is due to the emission from the jet lines. The Balmer 
lines originate at locations with almost no or only very low ve- 
locities. All of the observed flux is symmetrically distributed 
around the center of mass extending to velocities of roughly 
350 kms- 1 . 
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Fig. 10. Zoomed grey scaled representations of continuum 
subtracted, photometrically calibrated, high resolution spectra 
of H a , Hp, H 7 (left side), He n A 4686 and He n A 4542 (right 
side) in v, ^-coordinates. Not covered phases are marked with 
empty rows. In H Q the jets are clearly visible and show an or- 
bital motion similar to the Hell emission lines, whereas the 
main emission shows almost no orbital motion. The blue side 
of the Balmer lines are truncated by the P-Cygni absorption. 
In H 7 complex substructures can be detected. The fragmented 
line appearing in the H 7 spectra at a velocity of 1200 kms -1 is 
a remnant of a night sky line. Within the Helium lines a weak 
S-wave component is visible. 



The intense emission maximum is located near the center of 
mass. Part of this emission might come from the irradiated sec- 
ondary, but as the maximum of the emission is clearly shifted 
to Vy ~ kms -1 an irradiated secondary probably cannot ac- 
count for all of the observed flux. The Doppler map of Hp looks 
similar. Again, the whole flux including the intense maximum 
is distributed symmetrically around the center of mass. 



The Doppler maps from the Helium lines are slightly differ- 
ent. Comparable to the Balmer emission the observed flux orig- 
inates at locations with very low velocities. The intensity dis- 
tribution as a whole is distributed roughly symmetrical around 
the center of mass. But the emission maximum of the Helium 
lines is not centered at the center of mass. The coordinates 
of this intense maximum are shifted to approximately Vx ~ - 
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Fig. 11. Doppler-maps of H Q , Hp, He II A 4686 and He II A 4542 in velocity space (Vx, Vy), using the Fourier-filtered back- 
projection technique. The emission mainly originates in a "disk" centered at the center of mass, which extends to velocities of 
only 350 kms -1 . The maximum of the emission of the Helium lines is shifted to Vx ~ -100 kms -1 and Vy ~ -100 kms -1 . 
The dim ring in H a is caused by the jet emission. From the distribution of the emission in the Doppler maps it is clear, that the 
emitting material is not located in an accretion disk. 



100 kms -1 and Vy « -100 kms -1 in both maps. This emission 
maximum does not coincide with the gas stream trajectory. But 
due to its position in velocity space the emission might origi- 
nate in the elevated and radially extended accretion disk rim as 
proposed by Meyer-Hofmeister et al. (1998). Thus, the Helium 
emission might partly be related to the outermost accretion disk 
rim. The Hell (A 4686) map also indicates emission from re- 



gions near the white dwarf, which supports our assumptions in 
Sect. 3.3 for the estimation of the mass of the secondary. 

Another indication about where this material might be lo- 
cated is given by the fact that the emission is roughly centered 
on the center of mass and not on the primary as would be ex- 
pected for a typical accretion disk. This material might perhaps 
orbit the binary system around the mass center far outside of 
both Roche lobes. In a very simplified estimate one finds a dis- 
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Fig. 12. Kepler velocity at the rim of the accretion disk at 
0.8Rl (Rl is the volume radius of the Roche lobe of the white 
dwarf) over mass ratio q for various masses Mwd of the white 
dwarf and inclinations i = 50° and 90°. For the observed ve- 
locities it is clear that the emission lines cannot originate within 
the accretion disk. For medium inclinations an origin only at 
the rim of the disk is possible. 

tance of 10 119 cm for the above given binary parameters and a 
supposed velocity of 150 kms -1 for the emitting circum-binary 
material. This is about 3.5 times the distance of the binary sep- 
aration. 

4. Discussion and summary 

We performed high-speed photometric and spectrophotomet- 
ry observations as well as high-resolution spectroscopy of the 
bright galactic SSS RXJ0019 between 1992 and 1997. Our de- 
tailed optical studies result in a new determination of the orbital 
period. We observed a color variation in R — I symmetrical 
to the main minimum with an amplitude of Am = 0.15 m . The 
observed lightcurves change between having a flat appearance 
(quiet optical state) and showing humps and steps with a 1.8 h 
period (excited optical state). RXJ0019 can change between the 
two states from night to night. These variations might be due 
to short-term changes of the mass accretion rate from the com- 
panion star. 

Our spectroscopic investigation shows the well known 
emission line spectrum of mainly Balmer and Helium II lines. 
Symmetric emission lines adjacent to the Balmer lines show 
the presence of high velocity outflows (jets) probably originat- 
ing near the white dwarf. The jet lines show an orbital Doppler 
motion comparable to that of the He II lines. The velocity of 
the jet is quite low, indicating a medium to high inclination of 
the system. 

From the radial velocities we calculated the mass function. 
For the assumed range of medium to high inclinations we de- 
rive a low mass for the secondary star (0.3 < M 2 < 0.5Mq). 
Low mass companion stars might be common among the SSS. 
This might also be an explanation for the absence of any spec- 



tral features of the companion star in the high resolution spec- 
tra. 

We also observe P-Cygni profiles in the Balmer lines show- 
ing an orbital modulation. The P-Cygni absorption almost 
disappears between cf) — Q.l to 0.9. Higher transitions of the 
Balmer series up to H i2 are only detected by their correspond- 
ing P-Cygni absorption. The velocity of the directed wind re- 
sponsible for the P-Cygni profiles is « 590 kms -1 with the ab- 
sorption wings extending to 900 kms -1 . 

Our trailed spectra show substructures of different compo- 
nents. As these components are very weak or were observed 
in lines with bad S/N ratio the locations of these components 
cannot clearly be located by the means of Doppler tomography. 
The shifted intensity maximum in the Helium maps at Vx ~ - 
100 kms -1 and Vy ~ -100 kms -1 might be due to a radially 
extended and elevated accretion disk rim, consistent with the- 
oretical models of SSS. But we clearly see that most of the 
line-emitting material is not located within the accretion disk. 
Due to the high accretion rate and viscous processes involved 
we assume that, almost the whole accretion disk is optically 
thick. Therefore, no emission from Balmer and Helium lines 
originating at the inner parts of the disk can be observed. 

We only see line-emitting material with very low veloci- 
ties. The emission distribution is mainly symmetrical around 
the center of mass. As material cannot be stationary at the cen- 
ter of mass the emission must originate at regions with very 
low velocities. Low velocities can be found far outside of the 
binary system. Therefore, we propose that this SSS has a cir- 
cumbinary cocoon or disk of hydrogen and helium responsible 
for the observed emission lines. 
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